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Abstract. We consider solutions of the focusing cubic and quintic Gross- 
Pitaevskii (GP) hierarchies. We identify an observable corresponding to the 
average energy per particle, and we prove that it is a conserved quantity. We 
prove that all solutions to the focusing GP hierarchy at the L^-critical or L^- 
supercritical level blow up in finite time if the energy per particle in the initial 
condition is negative. Our results do not assume any factorization of the initial 
data. 



1. Introduction 

The mathematical analysis of interacting Bose gases is a rich and fascinating 
research topic that is experiencing remarkable progress in recent years. One of 
the fundamental questions in this field concerns the mathematically rigorous proof 
of Bose-Einstein condensation; for some recent landmark results in this direction, 
we refer to the works of Lieb, Seiringer, Yngvason, and their collaborators which 
have initiated much of the current interest in the field, see [21 EQI EH [22] and the 
references therein. 

Another main line of research focuses on the effective mean field dynamics of 
interacting Bose gases. In the recent years, remarkable progress has been achieved 
in the mathematically rigorous derivation of the nonlinear Schrodinger (NLS) and 
nonlinear Hartree (NLH) equations as the mean field limits of interacting Bose 
gases. For some recent fundamental results in this area, we refer to the works of 
Erdos, Schlein and Yau in [51[Tni[II], and also [Tni[IHl[25] and the references therein; 
see also [D [3 [51 [H [H [H [HI [M [13 [23 • 

The strategy of [Hllinilll] involves the following main steps, which are presented 
in more detail below: Based on the Schrodinger evolution of the given iV-body 
system of bosons, one derives the associated BBGKY hierarchy of marginal density 
matrices. Subsequently, one takes the limit N oo, whereupon the BBGKY 
hierarchy tends to an infinite hierarchy of marginal density matrices referred to as 
the Gross-Pitaevskii (GP) hierarchy. Finally, one proves that for factorizcd initial 
conditions, the solutions of the GP hierarchy are factorizcd and unique, and that 
the individual factors satisfy the NLS or NLH, depending on the definition of the 
original A^-body system. In the work at hand, we will focus on the case linked to 
the NLS. 

It is well known that for focusing L^-critical and L^-supercritical NLS, nega- 
tivity of the conserved energy implies blowup of solutions in H^. This is usually 
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proven by use of energy conservation combined with a virial identity, a method 
often referred to as Glassey's argument. In this paper, we are especially interested 
in the phenomenon of blowup of solutions for the GP hierarchy without assuming 
factorization of the initial conditions. More precisely, here we obtain an analogue 
of Glassey's argument for the GP hierarchy, and thereby, we establish blowup of 
solutions to the GP hiearchy under the condition that the initial energy is negative. 

First, for the convenience of the reader, we outline below the main steps along 
which the defocusing cubic NLS is derived as the mean field limit for a gas of bosons 
with repelling pair interactions, following [9l[ini|Tl]. For repelling three body in- 
teractions leading to the defocusing quintic NLS, we refer to [5]. We remark that 
it is currently not known how to obtain analogous results for the case of attractive 
interactions. 



(i) From N-body Schrddinger to BBGKY. Let 7/1^ G L'^[W^^) denote the wave 
function describing N bosons in M'^. To account for the Bose-Einstein statistics, 
i/iAT is invariant with respect to permutations tt G Sn, which act by interchanging 
the particle variables, 

i'N{x-^{l),X^(2), •■•,2^7r(Af)) = i>N(xi,X2, ■■■,Xn) ■ (1-1) 

We denote Ll{R'^^) -.^ {ip^ e L^{R'^^)\^pN satisfies (HH)}. The dynamics of the 
system is determined by the Schrodinger equation 

idtiJN = HniIjn ■ (1-2) 

The Hamiltonian Hff is assumed to be a self-adjoint operator acting on the Hilbert 
space LliR'''^), of the form 

^ 1 

77Ar = ^(-A,J + - VNix,~x,), (1.3) 

j=l l<i<j<N 

where Vn{x) = N'^^V(N'^x) with V E W^'''{R'^) spherically symmetric, for some 
suitable r, s, and for (3 G (0, 1) sufficiently small. 

The limit — s- oo is obtained in the following manner. One introduces the 
density matrix 

7Ar(t,X^,X^) = 1pN{t,Xj^)4'N{t,x'j^) 

where Xj^j — (xi, a;2, xat) and — {x[,x'2, .-.jx'pj). Moreover, one introduces the 
associated sequence of fc-particle marginal density matrices 7^'' (t) , for k — 1, . . . , N, 
as the partial trace of 7Ar over the degrees of freedom of the last {N — k) particles, 

7^' = Trfe+i,fc+2,...,A'|'0Af)(V'Ar| • 
Here, Trfe+i fc-|_2,...,Af denotes the partial trace with respect to the particles indexed 
by fc -I- 1, /c + 2, TV. Accordingly, is defined as the non- negative trace class 
operator on Ll{R'^^) with kernel given by 

iN^SLkillk) = j dxj^_i.^N{xj^,Xj^-ki3lkT3LN-k) 

= j dx^_ktpN{xk,Xj^^k)i'N{x^k^x^-k) ■ (1-4) 
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It is clear from the definitions given above that j'^'^ — Trfe+i7^^^'', and that 
Tr7^'' = ||V'w|li2(R.«) = 1 for aU iV, and all fc = 1, 2, N. 

The time evolution of the density matrix is determined by the Heisenberg 
equation 

%daN{t)^\HN.lN(t)\. (1.5) 

which is equivalent to 
idaN{t,xj^,x'j^) = -{A^ - A^' )-fN{t,Xj^,x'j^) (1.6) 



T7 i^N{xi~Xj)-VN{x'^-~Xj)]'jN{t,x^,x!j^), 



N 

l<i<j<N 

expressed in terms of the associated integral kernel. Accordingly, the /c-particle 
marginals satisfy the BBGKY hierarchy 

ida^''\t,x,;x:,) = -(A^^_-A^,)7W(i, 01^,4) 

+ ^ E [VNix,-x,)~VNix[^xr)h(''\t,x,;x:,) (1.7) 

l<i<j<fc 

N — k , ^ f 

H Jj—y^, / dxk+i[VN{xi - Xk+i) - VNix'i - Xk+i)] (1.8) 

J''''^^Ht,x^^,Xk+l■,x^^,x'k+l) 
where A^,^ := X]j=i ^2:^1 and similarly for A^'^. We note that the number of terms 
in p.7p is « ^ — > 0, and the number of terms in (jl.Sp is ''^^^'^'i as — > cx). 
Accordingly, for fixed fc, (jl.7p disappears in the limit N 00 described below, 
while (11.811 survives. 



(ii) From BBGKY to GP. It is proven in [9l[T0l[TT] that, for a suitable topology on 
the space of marginal density matrices, and as A^ — > cxd, one can extract convergent 
subsequences 7^'' — > 7'*^' for k E N, which satisfy the infinite limiting hierarchy 

ida^''\t,Xk;x:,) = -(A,^, -A,,)7«(t,a:fe;xy (1.9) 



k 

■6o^Sj-fc+i7'=+i(t,Xfe;x'fc), 



which is referred to as the Gross- Pitaevskii (GP) hierarchy. Here, 

(BM+n*+')(i.ais;) 

'<i««<i4.,l*fe -4.,)- 

7^''+^Ht,^fe,a;/c+i;x'fe,4+i), 

and bo — J V{x)dx. The interaction term here is obtained from the limit of (|1.8 
as A^ ^ 00, using that VNix) bQ5{x) weakly. We will set 60 = 1 in the sequel. 



(Hi) NLS and factorized solutions of GP. The link between the original bosonic A^- 
body system and solutions of the NLS is established as follows. Given factorized 
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/c— particle marginals 

k 

at initial time t = 0, with (po G H^{M.'^), one can easily verify that the solution of 
the GP hierarchy remains factorized for alH G / C M, 

k 

if G H^iR'') solves the defocusing cubic NLS, 

idt^ = -A,^ + \<l>\^(l>, (1.10) 

for t G /, and (/)(0) = 0o e 

Solutions of the GP hierarchy are studied in spaces of fc-particle marginals with 
norms IhWllL := Tr(S'W7W) < oo or he^^ll^i := (Tr(5('=)7W)2)i/2 < oo ^here 

^.(/c) _ nj=i(Va;,)(V:r;.), and i?^ = i/"(M'"= x M'^'^) for brevity. While the exis- 
tence of factorized solutions can be easily obtained, as outlined above, the question 
remains whether solutions of the GP hierarchy are also unique. 

The proof of uniqueness of solutions of the GP hierarchy is the most difficult 
part in the program outlined above, and it was originally accomplished by Erdos, 
Schlein and Yau in UHl |TT] by use of sophisticated Feynman graph expansion 
methods. In Klainerman and Machedon proposed an alternative method for 
proving uniqueness based on use of space-time bounds on the density matrices and 
introduction of an elegant "board game" argument whose purpose is to organize 
the relevant combinatorics related to expressing solutions of the GP hierarchy using 
iterated Duhamel formulas. For the approach developed in |19j . the authors assume 
that the a priori space-time bound 

\\B,,,+a^''+'^\\LiHl<C\ (1.11) 

holds, with C independent of k. The authors of [T^ proved that the latter is indeed 
satisfied for the cubic case in d = 2, based on energy conservation. 

Non-factorized solutions of focusing and defocusing GP hierarchies. As mentioned 
above, it is currently only known how to obtain a GP hierarchy from the N oo 
limit of a BBGKY hierarchy with repulsive interactions, but not for attractive 
interactions. 

However, in the work at hand, we will, similarly as in f7^, start directly from the 
level of the GP hierarchy, and allow ourselves to also discuss attractive interactions. 
Accordingly, we will refer to the corresponding GP hierarchies as cubic, quintic, 
focusing, or defocusing GP hierarchies, depending on the type of the NLS governing 
the solutions obtained from factorized initial conditions. 

Recently, in [7] , two of us analyzed the Cauchy problem for the cubic and quintic 
GP hierarchy in R'^, d > 1 with focusing and defocusing interactions, and proved 
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the existence and uniqueness of solutions to the GP hierarchy that satisfy the space- 
time bound which was assumed in . As a key ingredient of the arguments 
in [7] a suitable topology is introduced on the space of sequences of marginal density 
matrices, 



© = {r = (7W(a;i,...,a;fc;a;'i,...,4))fceN|Tr7W < oo } . (1.12) 



parametrized by ^ > 0, and the spaces = {T G (S|||r||-Ha < oo} were intro- 
duced. The parameter ^ > is determined by the initial condition, and it sets the 
energy scale of a given Cauchy problem; if F g H'^, then is an upper bound 
on the typical ff"-energy per particle. The parameter a determines the regularity 
of the solution. In [7j, the local in time existence and uniqueness of solutions is 
established for cubic, quintic, focusing and defocusing GP hierarchies in Ti.'^ for a in 

a range depending on d, which satisfy a spacetime bound ||Sr||Ljgj--H° < C'||ro||-H° 
for some < ^ < (here 13T { B^;_^_p"f^''^i'^ )keN)- The precise statement and 
the associated consequences that we will use in this paper are presented in the next 
section. This result implies, in particular, (|l.lip . 

In this paper we study solutions of focusing GP hierarchies without any factor- 
ization condition, and especially establish the following results characterizing the 
blowup of solutions: 

(1) For defocusing cubic GP hierarchies in d — 1,2,3, and defocusing quin- 
tic GP hierarchies in d = 1,2, which are obtained as limits of BBGKY 
hierarchies as outlined above, it is possible to derive a priori bounds on 
llT^'^HOIIij^/fi based on energy conservation in the A^-particle Schrodinger 
system, see [SJ [13 E], and also [SJ [H]. However, on the level of the GP 
hierarchy, no conserved energy functional has so far been known. We iden- 
tify an observable corresponding to the average energy per particle, and we 
prove that it is conserved. 

(2) Furthermore, we prove the virial identity on the level of the GP hierar- 
chy that enables us to obtain an analogue of Glassey's argument from the 
analysis of focusing NLS equations. As a consequence, we prove that all 
solutions to the focusing GP hierarchy at the L^-critical or L^-supercritical 
level blow up in finite time if the energy per particle in the initial condition 
is negative. 

Organization of the paper. In Section [5] we present the notation and the pre- 
liminaries. The main results of the paper are stated in Section [31 In Section [¥| we 
identify the average energy per particle and prove that it is a conserved quantity. 
In Section [SI we derive a virial identity that enables us to prove an analogue of 
Glassey's blow-up argument familiar from the analysis of NLS. The analogue of 
Glassey's blowup argument is presented in Section [S] 



It is determined by the generalized Sobolev norms 




(1.13) 
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2. Definition of the model and preliminaries 
We introduce the space 

oo 

© := 0L2(R'''= X R'^'') (2.1) 

k=l 

of sequences of density matrices 

r:-(7('^WN (2.2) 

where 7*^*^^ > 0, Tr7*^'^^ — 1, and where every J^'^^X)., x'lJ is symmetric in all 
components of Xj,, and in all components of xj^, respectively, i.e. 

7('')(a;^(i), ...,x^(fc);<,(i), ...,a;^,(fc)) = i'^''\xu ■■■,Xk;x[, (2-3) 
holds for all tt,7t' £ Sk- 

Moreover, the fc-particle marginals are hermitean. 



7^'^fefc;2i;) = 7«(2^;^fe)- (2.4) 

We call r = (7(''))a;gn admissible if 7^^) = Tri,+i,._fc+p7(''+5), that is, 
7('H£.;2'fe) (2.5) 
dxfc+i . --dXk+R 7(''+^'(xfc,a;fe+i, . . . , x^+e ; x'j., x^+i, . ..,Xk+^) 
for all fc e N. 

We will use the following convention for the Fourier transform, 
l{xk:ilk)= J rflifc^^4e-'=-'^--'^-'^7(^fc;u'fe). 

Let < 5 < 1. We define 

:= [r e &\\\r\\n^ < ^} (2.6) 

where 

00 

l|r||«^" = E^'ll^^'^ll«°(K'"=xK-)' (2.7) 

fe=l 

with 

II 7^''Ml//°(R'»=xK<"') = II 5'^'''"^^'''' llL2(R<ifcxR'"=) : (2.8) 

and S'^'^'") := n^=i(^2:j)"(^2:' Clearly, H'^ is a Banach space. Similar spaces 
are used in the isospectral renormalization group analysis of spectral problems in 
quantum field theory, [4J. 

Next, we define the cubic, quintic, focusing, and defocusing GP hierarchies. Let 
p e {2,4}. The p-GP (Gross-Pitaevskii) hierarchy is given by 

k 

ida^"^^ ^ ^[-A,^,7«] + /iBfe+|7^'+^^ (2.9) 
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in d dimensions, for fc G N. Here, 

(i?fc+|7^'+*^) (t,a^i,...,Xfe;x;,...,4) (2.10) 

k 

■= 51 (-Si;fc+i,...,fc+£7^''^^^) (i,a;i,...,Xfe;xi,...,4) 

A; 

lit (^l;fc+i,...,fe+f7^''+^^) (^,2;i,...,a;fc;x'i,...,x'fe) 

where 

(^l;fc+i....,fc+f 7^*-'+*^) {t,xi,...,xk;x[,..., 4) 
fc+f 

-a;£)(5(xj - a;;,)7(''+5)(t, xi, . . . , x^+i ; x'l, . . . , 4+e)> 

and 

(^'fc+i,....fc+f 7^''+*^) (t, . . . , Xfc; x'l, . . . , 4) 

= J dxk+i---dxk+Edx'f,j^i---dx'^^E 
fc+f 

J]^ -Xi)S{x'^ -4)7(''+^)(i,a;i,...,Xfe+£;a;'i,...,4+£). 

The operator i3fc_|_£7''^''"5) accounts for | + 1-body interactions between the Bose 
particles. We note that for factorized solutions, the corresponding 1-particle wave 
function satisfies the p-NLS idtcj) ~ —Acf) + fi\(f>\P(j). 

We refer to (|2.9p as the cubic GP hierarchy Hp — 2, and as the quintic GP 
hierarchy if p = 4. Also we denote the L^-critical exponent by p]^2 — ^ and refer to 
(j2.9p as a -critical GP hierarchy if p = p]^2 and as a -supercritical GP hierarchy 
if p > pl2. Moreover, for /i = 1 or /i = — 1 we refer to the GP hierarchies as being 
defocusing or focusing, respectively. 

To obtain the blow-up property, we need a result providing a blow-up alternative. 
This is a usually obtained as a byproduct of the local theory. In the context of the 
local theory developed in [7], we recall the following two theorems: Theorem 12.11 
which establishes the local well-posedness of the GP equation and Theorem [^3] that 
gives lower bounds on the blow-up rate. In order to state these two theorems, we 
recall that in [7] the GP hieararchy was rewritten in the following way: 

idtT + A±r = ^i?r, (2.11) 

where 

A±r := ( Ai'^^W )feeN with A^"^ = A,^ - A,;, 
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and 

5r (i?,+ z7('=+i)),eN. (2.12) 

Also the following set 2l(c!,p) was introduced in f7], for p = 2, 4 and d > 1, 

( (5,00) ifd=l 
nd,p) = I (I -21^'^) ifd>2and(d,p)^(3,2) (2.13) 
[ [1,00) if (d,p) = (3,2) 

which we wih use to account for the regularity of solutions. 

Now we recall the local well-posedness theorem proved in [7| . 

Theorem 2.1. Let ^1 > 0. Assume that a £ 2l((i,p) where d>l and p G {2,4}, 
and < J] < I sufficiently small. Then, the following hold. 

(i) For every Fq G T^I*^ , there exist constants T > and < ^2 ^1 such that 
the following holds. There exists a unique solution T{t) in the .space 

{r e ||Br|Lj^^„„^ < 00} . 

where I := [0,T]. In particular this solution satisfies the Strichartz-type 
bound 

l|sr|Uj^^«»^ < c(r,d,p,a,6)l|ro||«^^. (2.14) 



(ii) The uniqueness of solutions in the space L'^jTi.'^^ is characterized as follows. 
Given Tq G ^j^; assume that there are constants T > and < ^2 < 
such that there exists a solution T{t) of the p-GP hierarchy Ii2. in the 
space L'^jH^^ . 

Then, this solution is unique in L'^jTi.^^ if and only if || BT < 00 

holds for some ^ > 0. 

Definition 2.2. We say that a solution T{t) of the GP hierarchy blows up in finite 
time with respect to H°' if there exists T* < 00 such that the following holds: For 
every ^ > there exists T^y ^ such that ||r(t)||-^Q 00 as t T^t- Moreover, 
T*j. / T* as C ^ 0. ' 



For the study of blowup solutions, it is convenient to introduce the following 
quantity. 



Definition 2.3. We refer to 

AvH»(r) [sup{e>0|||r||«o <oo} 



(2.15) 



AvL.(r) sup{e>o| ||r|U.. < ^} 



(2.16) 



respectively, as the typical (or average) H°' -energy and the typical U" -norm per 
particle. 
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We then have the foUowing characterization of blowup. 

Lemma 2.4. Blowup in finite time of T{t) with respect to as t /' T*, as 

characterized in Definition \2.2l is equivalent to the statement that Av h" {^{t)) ^ oo 
as t y T* (and similarly for ). 



Proof. Clearly, AvHo(r) is the reciprocal of the convergence radius of ||r||>^^ as a 
power series in Accordingly, ||r||->^o < oo for ^ < Av^/c (F)^-^, and ||r||-Hc. = oo 
for C > AvH- (r)-i. 



Blowup in finite time of T{t) in H°' as t /' T*, as characterized in Definition 



\2.2\ is equivalent to the statement that the convergence radius of ||r(t)||-Ha, as a 
power series in ^, tends to zero as t y T*. Thus, in turn, blowup in finite time of 
r{t) with respect to H" is equivalent to the statement that Av^fQ(r(t)) — > oo as 
t/T*. □ 

The following theorem from [7] gives lower bounds on the blow-up rate. 

Theorem 2.5. Assume that T{t) is a solution of the (cubic p — 2 or p — A quintic) 
p-GP hierarchy with initial condition r(io) = Fq G H.'^ , for some ^ > 0, which 
blows up in finite time. Then, the following lower bounds on the blowup rate hold: 

(a) Assume that ^ < p < -j^^- Then, 

Thus specifically, for the cubic GP hierarchy in d — 2, and for the quintic 
GP hierarchy in d — 1, 

(Av^.(r(i)))^ > -^-71, (2.18) 

|<-r*|2 

with respect to the Sobolev spaces H°', H'^. 

(Avi.(r(t)))^ > /or^<r. (2.19) 



Remark 2.6. We note that in the factorized case, the above lower bounds on the 
blow-up rate coincide with the known lower bounds on the blow-up rate for solutions 
to the NLS (see, for example, [51). 

We note that 

r - (|0)(</)|«'=)feeN ^ AvH=(r) - UWj,^ andAvi.(r) = UWl^ (2.20) 
in the factorized case. 

The fact that T e means that the typical energy per particle is bounded by 
Av^fo (r) < Therefore, the parameter ^ determines the i/"-energy scale in the 
problem. While solutions with a bounded i?"-energy remain in the same Ti,'^ for 
some sufficiently small ^ > 0, blowup solutions undergo transitions 7^1^^ TL'^^ 
~^ " ' where the sequence £,i> £,2> ■ ■ ■ converges to zero ast^T*. 
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We emphasize again that (AvAr(r)) ^ is the convergence radius of ||r||A/'j as a 
power series in ^, for the norms N = L"^ and A/"^ = T-C^, C^, respectively. 

3. Statement of the main results 

The foUowing two Theorems are the main results of this paper. First we prove 
energy conservation per particle for solutions T(t) of the p-GP hierarchy. More 
precisely, 

Theorem 3.1. Let < ^ < 1. Assume that r{t) e H'^ , with a > 1, is a solution of 
the focusing (fi = — or defocusing (fi — I) p-GP hierarchy with initial condition 
Tq G Ti'^ . Then, the following hold. Let 

k 

Ekim) iTr(^(-A.j7"=nO) (3.1) 

i=i 

+ ^^(E4..i,...,..f7<^"^'w). 

Then, the quantity 

Ene(rW) := Y.^'Mm) (3.2) 

k>l 

is conserved, and in particular, 

En^im) = {J2ke)E,iT{t)), (3.3) 

k>l 

where X]fe>i ^C'^ °° '^^V < ^ < 1. 

We recall that Tr(^) means integration of the kernel A{x, x') against the measure 
/ dxdx' 6{x ~ x'). We note that for factorized states T{t) = {\(f>{t)){(j){t)\'^*')ke¥ii one 
finds 

Eiim) = ^l|V0(t)||i. + ^||<^(t)||^+^., (3.4) 

which is the usual expression for the conserved energy for solutions of the NLS 
i9t0 + A(/) + ^|0|P(/) = O. 

Theorem 3.2. Let p > pl^. Assume that T(t) — {^^'^^t) )keN solves the focusing 
(i.e., /i = —1) p-GP hierarchy with initial condition T(0) G 7i| for some < ^ < 1, 
with Tr( 2:27(1) (0)) < 00. If Ei{T{0)) < 0, then there exists T* < 00 such that 
Avjji(r(t)) -^00 ast /T*. 

Remark 3.3. We note that Theorem is proved under the assumption that 
Tt{x^"/'--^"' (0) ) < 00, which is analogous to the finite variance assumption in the 
case of Glassey's blow-up argument for the NLS (see, e.g. [5] Theorem 6.5.4). 

As a motivation for the proofs presented below, we briefly recall the application 
of Glassey's argument in the case of an L^-critical or supercritical focusing NLS. 
We consider a solution of idtcj) = —A(j} — with </)(0) = 4>o € _ffi(R'') and 

P > Pl^ = 2' ^^^^ tli^t the conserved energy satisfies E[(j){t)] := 5II V(/)(t)|||2 — 
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^^II<?^(^)IIlp+2 = E[(f>Q] < 0. Moreover, we assume that || |a;|0o < oo- Then, one 
considers the quantity V(t) := {4>{t),x'^(f){t)), which is shown to satisfy the virial 
identity 

d^Vit) = 16i?[0o] - idP^^Wmiltl^ ■ (3.5) 

Hence, if -E[0o] < 0, and p > Pl^, there exists a finite time T* such that the positive 
quantity V{t) \ as t y T* . Accordingly, this imphes that ||0(i)llHi(R'*) oo as 
t T* (for more details, see Section[6]). This phenomenon is referred to as negative 
energy blowup in finite time for the NLS. In the sequel, we will prove analogues of 
these arguments for the GP hierarchy. 



4. Conservation of energy 

In this section, we prove Theorem 13.11 We first demonstrate the proof for the 
cubic case, p — 2. To begin with, we note that 

Ek{T{t)) = kE^{T{t)) (4.1) 

for V{t) = (7(^-)(0)fc gN a solution of the cubic GP hierarchy, where by definition, 
all 7''^)'s are admissible. To prove this, we note that (|3.ip can be written as 

k 



Eu{T{t)) = ^[lTr((-A.j7('H0) 



+ ^Tr(i?i,+i7^'+''(i))], (4.2) 

where each of the terms in the sum equals the one obtained for j = 1, by symmetry 
of 7*^*^^ and 7('^+i) with respect to their variables. We present the detailed calcula- 
tion for the interaction term, and note that the calculation for the kinetic energy 
term is similar. Consider 1 < i < j < k. We have that 

(4fe+i^^'^'')(^i' =^2, Xfe; x'l, 4, 4) (4.3) 

' ^ {x\ , X2 1 .. ., , . .. , Xk , Xj , Xi , X2 , . .. , , . .. , Xj- , .. ., Xf^ ■} '^j) J 

and 

iBl,+,j^''+'^){xuX2,...,Xk;x[,x'„...,xi) (4.4) 

— 'J^- ' ^ {xi , . .. , , . .. , , .. Xk 1 Xif Xi , . .. , X^ , . .. , X^- , . .. , X^ J Xi') . 

Thus 

Tr(4,+i7('+')) (4.5) 

7^ ^ (^1 ; . . . 5 Xi , . . . , Xj , . . . , , Xj Xj , . . . , , Xj^clX]^dX2 • • - dxf^ . 



By the symmetry of 7^'^^^-'(2:j,_|_]^; a;'fc+i) with respect to the components of 
and aijj+i, respectively, 

Tr(i?l..+i7('+'^)-Tr(i?i,^.i7('=+^)) 
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for all «, j e {1, . . . , k}. Thus, 

Mm) = fc[iTr((-A.j7<'Ht)) 

+ ^TTiBl,^,^(^+'Ht))] (4.6) 

follows. 

We can go one step further and note that 

(^i,fe+i7*''^^')(a;i,a:2, Xfe; x'l, 4, ■■■,x'k) 

= 'y^^^^\xi,X2,-;Xk,xi;x[,x'2,...,x',„xi). (4.7) 

Then 

^'•''^^^xi,xi,X2, ...,Xk;xi,xi,X2, ...,Xk)dxidx2-.-dxk 

where in the last equality we used symmetry of 7^'^+^) . On the other hand, 

(b;,T™)(ii;i',) = ^">^Xl,Xl■I[,Il) (4.8) 

■..'«)(.„.„.,....,..;.;..„.„...,..)<ix....<i.. 

by repeated use of the admissibility of 7'^^-', for j = 2, . . . , k + 1. Thus, 
and 

Tr(Sj,,+i7('=+^^) = Tr(i3,V^)) 

follows. Therefore, 

i?A^(r(i)) = k[^Trii-A,,)j('Ht)) + ^TTiBl,^('\t))_ 

= fci?i(r(t)), (4.9) 

as claimed. 

The fact that p.3p then follows is evident. 

Next, we verify that Ei ( r(i) ) is a conserved quantity, which means that Ei ( r(t) ) 
Ei{T{0)) for aU t e R. 

For the proof, the following auxiliary identities are very useful: 
dxdx'S{x — x')Vx ■ Vj;' A(x; x') 

dxdx5{x - x')/^xA{x;x') (4.10) 

dxdx'5{x - x')l^^'A{x]x') . (4.11) 
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To prove ()4.10p . we note that 

J dxdx'S{x — x')Vx ■ \/x'A{x; x') 

= - j dudu' j dxdx'6{x - x') e*"^'*"'^' u ■ u'A{u; u') 
= — J dudu'S{u — u')u ■ u'A(u;u') 
= — j dudu'diu — u') V? A{u] u') 

dudu' j dxdx'5{x - x') e™^-™'^' u^A{u; u') 
dxdx5{x~x')l:^,j,A{x;x'). (4.12) 



The equahty (|4.1ip can be proved in a similar way. 

We now return to the proof of Ei{V{t) ) = i?i(r(0) ). For fc = 1, we consider 

7*^"'")(x, x') where 

Ei{r{t) ) = -iTr(V, -V^-tW) (4.13) 

in symmetrized form. Here, we have introduced the shorthand notation 

S{xi = X2 = x[ = x'2) '■= i{x\ — x-2)^(x\ — x'^8(xi — x'2) ■ (4-14) 

Clearly, 

idtEiiTit)) = (/) + (//) + (///) + {IV) (4.15) 

where 

(/) := iTr( V, • V.KAx - A,07^'^ ) • (4.16) 

(//) -|Tr(V,-V,. (51,27^2))) (^^^7) 

(///) := — ^ J dxidx2dx'idx'2 S{xi = X2 = x'l = X2) (4-18) 

(A^i + A^^ - Aj.'^ - A^'Jj^^^xi, X2; x[, x'2) 

{IV) :— ^ J dxidx2dx'idx'2 S{xi ^ X2 = x'^ — x'^} 

(B37^'^)(a;i )■ (4-19) 
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The term (/). We claim that (/) = 0. We note that 

= j dudu' j dxdx'6{x -x')u- ^'(w^ - (u'f) e™=^-™'^' jW (u; u') 
= j dudu' 5{u -u')u- u'{u^ - {u' f ) 7(1) {u; u') 

= . (4.20) 
This proves the claim. 

The term {IV). We claim that {IV) also vanishes. Indeed, 

{IV) = ^ J dxidx2dx'idx'2 6{xi = X2 = x'l = x'2) 

{B3'y^^^){xi,X2;x[,x'2) 

= Y.U^-^dX2dxW{x,=X2=x[=x',) 

xi) -'y^^\xi,X2,x[;x[,X2 x[) 
+ 7^^) (a;i , a;2 , a;2 ; x[ , ^2 , a;2 ) - 7(^^2:1 , 2:2 , ^2 ; x[ ,x'^,x'<^ 



= 0. 

This proves the claim. 



(4.21) 



The term {III) . By symmetry of 7^^^ {xi ) in {xi,X2), and in (xi,X2), we 

S^d 

{III) = 



-2 ^ y dxidx2dxidx2 6{xi = X2 = x'^ = x'2) 



(A^i - Aa;'j7(2)(a;i,a;2;a;i,a;2) 
= ~ j dxxdx'-^5{xi-x'^) (4.22) 

( A^, 7^^) ) (a;i , x'l ; a;'i , x'l ) - ( A^/ 7^^^ ) (a;i , xi ; a;i , ) 

We will show that this term is canceled by the term (//). 



The Icnii [II). We have 



{II) 



Tr(V,-V,-(Bi,27^2))) 



(4.23) 



= dxdx'5{x-x')'Wx-'^x'{l^'^\xi,xi;x'^,x-i)-^^'^^{xi,x'^;x'^,x'i)) 
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(//) = I dxidx[S{xi-x[) (4.24) 



Now we use (|4.10p . and we obtain 
2 . 

{^x[l^'^^){xi,xi;x[,xi) - {A^^'^^'^^){xi, x[; x[, x[)) 

Hence (IT^ and (PTM]) imply (//) = -(///). 

We conclude that dtEi{r{t)) 0. Therefore, Ei{T{t)) = Ei{T{0)) is a conserved 
quantity. It represents the average energy per particle. For the quintic case p = 4 
(or similarly in more general cases p G 2N) , the above arguments can be adapted 
straightforwardly. □ 

5. ViRIAL IDENTITIES 

In this section, we prove the virial identities necessary for the application of a 
generalized version of Glassey's argument. According to our previous discussion, it 
is sufficient to consider only 7'^^^(xi; x'^^). 

5.1. Density. In what follows, we drop the superscript "(fc)" from 7''''). It will be 
clear from the number of variables what the value of k (in this part of the discussion, 
fc = 1 or fc = 2) is in a given expression. 

We write 

7(x;a;') = J dwdw'e™^-™'^'7(w; w') (5.1) 

and define 

p{x) -/{x]x) = j dvdv' e'^''-''">''^{v;v'). (5.2) 

Thus, 

dtp{x) = j dvdv'e''-^-'''^'=da{v;v') (5.3) 
= -i J dvdv'e'^^-^">^ {A,'^^,)j{v,v') 
+ - f dvdv' e'^"-'"'^'-'1hXi{v;v'). 



First we notice that 

dvdv'e'^^-^'^-{Aj^,,)j{v,v') 
= ^J dvdv' e'("-"')"(v2 -(«')') 7(w;w') 

dvdv' e'('^-"')^(w + w')(^'-i^')7(«;«') 
= V^- I dvdv' e'^''-'"'y-'{v + v')^{v]v'). (5.4) 
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On the other hand, 

dudqdu'dg'e*("^+''^-"'^'-''y') 7(u, q; u' , q') 
dudqdu'dq'e^^^+'l-'i'^''-''''''^ 7(u, q; u' , q') . (5.5) 

Therefore, 

Bl^{v] v') ^ J dxdx' e-™^+*^'^' Bl^-/{x; x') 

dudqdu' dq' 6{u + q — q ~ v) 5{v' — u') 7(u, q; u' , q') 
^ dqdq'j{v-q + q',q;v',q'). (5.6) 

Likewise, one obtains 

Bl^{v;v') ^ j dxdx' Bl .^-i{x;x') 

dudqdu' dq' 5 [v — u) 5{v' — {u' + q' — q)) 7(u, g; u', q') 

^ dqdq'^{v,q;v' + q-q',q'). (5.7) 

Thus, 

i J dvdv' e^^"-^')^ S^^(w; w') 

. l/*..,V<-'NSr^,,./)-SS(.;.')) 

= i y" dvdv'dqdq' e^^"^"')^ 7(w - g + g', g; v' , q') 

- i y dvdv'dqdq' e'(^-"')=^ 7(w, g; v' + g - g', g') (5.8) 
= 0, 

where the last equahty is obtained by applying the change of variables v — > v — q + q' 
and u' — » u' — g + g' in the second term of (|5.8p so that the difference v — v' remains 
unchanged. 



Therefore, by combining (|5.3p . (|5.4[) and (|5.8|) we conclude 

5tp(x) - V, • P = 0, 

where 

corresponds to the momentum. 
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5.2. Morawetz action. We define 

M := j dxx-P . (5.10) 

The time derivative is given by 

dtM = J dxx-dtP - (Im) + (IIm), (5.11) 

where {Im) is the kinetic, and {IIm) the interaction term. 
We have 

{hi) = j dxx- j dudu' e'^''-'^'^'-'{u + u'){u^ -{u')'^)^{u;u') 

dudu'^{u;u') j dx x ■ [{u + u') ® {u + u')]{y ^e'^''-'''^'') 
<i„.„',(„;„')T.[(„ + „',«,„ + „')|/.....-0. 

du du' 5{u - u')Tr[(u + u') ®{u + u')] ^{u; u') 
= ^ j duu^^{u;u) . (5.12) 
which is 8 times the kinetic energy of one particle. 

5.3. Interaction term. Next, we study the interaction term 

{IIm) = ^i j dxx--^ j dvdv' e'''"""'^'' {v + v')B^{v]v'). (5.13) 

(A) The cubic case. 

As we have seen in (|5.5p . we have 

B;,.(...') = /^"We'<<— 
Therefore, by (|5.6p 

^i,27(w; v') = / dqdq' j{v - q + q\ q; v' , q') . 



Likewise by (15. 7p 

Bl2l{v;v')^ I dqdq'j{v,q;v' + q-q',q'). 
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Now we determine the term (IIm) in J dx x ■ dtP that involves the interaction. 
To this end, we first consider 

i Jdvdv' e^("-"')- {v + v')IhXf{v;v') 

= \ j dvdv' e*("-''')x ^ (BlPfiv; v') - Bl~^iv; v')) 
= \ j dvdv'dqdq' e*(''-"')^ {v + v') j{v ~ q + q\ q; v' , q') 

-- I dvdv'dqdq' e*("-''')^ + ^{^y^ q;v' + q- q' , q') . (5.14) 



In the last term, we apply the change of variables v — > v — q + q' and v' — > v' — q + q', 
so that the difference v — v' remains unchanged. We obtain that the above equals 

i J dvdv'dqdq' e*(^-"')^ {v + v')^{v - q + q' , q; v',q') 

- i y" dvdv'dqdq' (w + - 2g + 2g') 7(1; - g + g', g; w', g') 

= - J dvdv'dqdq' e'(''-"')^ 7(1; - g + 9', 9; v' , q') ( (u + v') - (w + u' - 2^ + 2g') ) 

= \ j dvdv'dqdq' e'(^'-"')^ 2(g - q') ^{v-q + q' , q; v' , q') . (5.15) 
The contribution of this term to the integral J dx x ■ dtP is given by 

— J dxx ■ J dvdv'dqdq' e''^^~^ 2{q ~ q')^{v ~ q + q' ,q;v' ,q') . (5.16) 
Next, we express everything in position space. 
We have that the last line equals 

- J dxx- J dXdYdX'dY' j dvdv'dqdq' e*("-^'):r 2{q - q') 

^i(-{v-q+q)X-qY+v' X' +q'Y') ( Y V X' Y'\ 



^l(-(v-q+q )JL-qY+v +q Y ) ry (^X ,Y ] X\Y 

^iv{x-X)-tv'(x~x') 2x.{q- q') e+^q(x-Y)-q'{X-Y') 



^ J dx J dXdYdX'dY' -f{X,Y;X',Y') J dvdv'dqdq 



= -fijdxj dXdYdX'dY'-f{X,Y;X',Y') J dqdq' 

S{x - X)5ix - X') 2X ■ Vx e+*«(^-^)-"?'(^-^') 

^ -I^J dXdYdY' -/{X,Y;X,Y') 

2X -VxSiX - Y)5{Y -Y') (5.17) 

= y" dXdY j{X, Y; X, Y) 2X ■ S{X - Y) 

= /i y dXdY 6iX - Y) {2d + 2X -Vx hiX, Y; X, Y) (5.18) 
where we have written 5{X - Y)6{X - Y') = 6{X - Y)6{Y - Y') to get (fSTT]) . 
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Now we note that 



J dXX ■ VxliX, X; X, X) (5.19) 
= J dXdYd{X -Y){X -S/x + Y ■S/Yh{X,Y;X,Y) 
= J dXdY 6iX -Y)iX ■ VxliX, Y;X,Y) + Y- VyliY, X; Y, X) ) 

= J dXdY6{X -Y){2X ■VxliX,Y;X,Y)) (5.20) 

where we used the symmetry ^{XjY; X,Y) — j(Y, X;Y, X), and renamed the 
variables in the last term. Clearly, (|5.19p equals 

-d J dX-f{X,X;X,X) (5.21) 

from integrating by parts. 

Therefore, combining (l5l^ . ([QO]) and (l53T|) 

{IIj^^j) = A* y dXdY S{X - Y) {2d + 2X -Vx h{X, Y; X, Y) 

= ^ y" dXdY 6{X - Y) {2d - d )7(X, Y; X, Y) 

= fid J dXj{X,X;X,X). (5.22) 
This is the desired result for the cubic case. 

(B) The quintic case. 



Now we give a sketch of the calculations related to the interaction term in the 
quintic {p — 4) case. Again it suffices to consider k = 1. 

Since in the case when p 4 we have 
5J,,37(.; .') = / dy' d^ d.'5{. -y)S{.- y') S{. - .) S{. - .') 

dudqdrdu' dq'dr'e^'''^+''y+''~''''''-'''y' ^{u, q, r; u\ q\ r') 

dudqdrdu' dq'dr'e'^^''+''+'-'''-''^''-'''''"> 7(u, q, r; u' , q' , r') , (5.23) 
by taking the Fourier transform we obtain 

BlZfriv.v') = [ dxdx'er^^^+^^'^' Bl^^,^{x;x') 

dudqdrdu' dq' dr' 6 {u + q + r — q' — r' ~ v) 

S{v' - u')j{u,q,r;u',q',r') 
dqdrdq' dr' ^{v ~ q — r + q + r' , q, r; v' , q' , r') . (5.24) 
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As in the cubic case, to determine the term (IIm) in / dxx-dtP, we first observe 
that 

i J dvdv' e^-"-"'y^ {v + v')B^-i{y;v') 

= i 1 dvdv' e^^-^'^ {v + v') (BQy(«; v') - Bf^^iv; v')) 
= i y dvdv' dqdq' drdr' e'^"'""'^'' {v + v')^{v - q ~ r + q' + r' , q, r; v' , q' , r') 
- - [ dvdv' dqdq' drdr' e^^""""')^ {v + v') j{v, q,r;v' + q + r - q' - r' , q' , r') , 



which after performing the change of variables v -^v — q — r + q' + r' and v' — > 
v' — q — r + q' + r' in the last term, becomes 

- J dvdv' dqdq' drdr' g^^""''')^ 2{q + r - q' - r') j{v - q - r + q' + r' ,q,r;v' ,q' , r') . 
Hence the contribution of this term to the integral J dx x ■ dtP is given by 

dxx- J dvdv' dqdq' drdr' e*'""" 2{q+r—q'—r') ^{v—q—r+q'+r, q, r; v' , q' , r') , 
which we express in the position space as follows 

- j dxx- j dXdYdZdX'dY'dZ' j dvdv' dqdq' drdr' e''^"'''"'^ 2{q + r - q' - r') 

^i{-{v~q-r+q'+r')X-qY-rZ+v'X'+q'Y'+r'Z')^^-^ y Z'X' Y' Z') 

= J j dx j dXdYdZdX'dY'dZ' -f{X,Y,Z-X',Y',Z') J dvdv' dqdq' drdr' 

2x-{q + r-q' - r') 

^+iqiX-Y)-iq'{X-Y') ^+ir(X-Z)-ir'(X-Z') 

= -^l j dx j dXdYdZdX'dY' dZ' -f{X,Y, Z; X' ,Y' , Z') J dqdq' drdr' 

Six - X)Six - X') 2X ■ Vx e+*«(^-^)+-(^-^)-^9'(x-i")— '(^-z') 

= -fJ-J dXdYdZdY'dZ' -/{X,Y,Z;X,Y',Z') 

2X ■ Vx S{X - Y) 6{Y - Y') S{X - Z) 6{Z - Z') (5.25) 

= -Ai y dXdYdZ -fix, Y, Z; X, Y, Z) 2X ■ Vx S{X - Y)6{X - Z) 

= ^ y" dXdYdZ 5{X - Y)5{X - Z) {2d + 2X -Vx )l{X, F ; X, Y) (5.26) 
where we have written 

5{X - Y)5{X - Y')5{X - Z)5{X - Z') = 5{X - Y)5{Y - Y')5{X - Z)5{Z - Z') 
to get ((05)) . 



BLOWUP SOLUTIONS FOR THE GP HIERARCHY 21 

On the other hand, usmg symmetry of 7(X, Y, Z; X, Y, Z) we obtain 
J dXX ■ Vxl{X, X, X; X, X, X) 

= J dXdYdZ 5{X - Y)5{X - Z) 

{X -Vx + Y -Vy + Z-yzhiX,Y,Z;X,Y,Z) 

= J dXdYdZ 5{X - Y)5{X - Z) 

( X ■ Vxl{X, Y, Z- X,Y,Z) + Y- Vyi{Y, X, Z; Y, X, Z) 

+ Z-Vzl{Z,Y,X;Z,Y,X)) 

= j dXdYdZ S{X - Y)6{X -Z){3X ■ Vx^fiX, Y, Z; X, Y, Z) .) (5.27) 
However, by the integration by parts, 

j dXX ■\/xj{X,X,X;X,X,X) = -d J dX j{X,X,X;X,X,X), (5.28) 
so by combining (|5.27p and (|5.28p we obtain 

j dXdYdZ S{X - Y)b(X ~Z){ X- VxliX, Y, Z; X, Y, Z) 

j dX-f{X,X,X;X,X,X). (5.29) 



d 

Therefore 



{11 m) 



/i J dXdYdZ 6{X - Y)S{X - Z) {2d + 2X-Vx )l{X, Y, Z; X, Y, Z) 
fi I dXdYdZ S{X - Y)5{X - Z) {2d - 2d/3) )-f{X, Y; X, Y) 

J dX 7(X, X,X;X, X, X) . (5.30) 



id 



3 

This is the desired result for the quintic case. 
( C) The general case p € 2N 



The above calculation on the interaction term can be reproduced for a general 
even number p and it that case we obtain: 

{Ihi) = dXidX2...dXi+E S{Xi - X2)...S{Xi - X,+e) 

{2d + 2Xi-Vx,h{Xi,...,X,+ .;Xi,...,X,+ p) 

= y dXidX2...dXi+p S{Xi ~ X2)...6{Xi - Xi+e) 

{2d - 2-^ )7(Xi, X1+5) 
^y'rfX7(X_X;X_£). (5.31) 
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Now we combine (|5.12p and (|5.3ip to conclude that: 

j dx x'^'-f{x, x) 

= 2 J dxx-dtP 

^ ^ j duu^^{u\u) + Ai^^ j rf^7(£^^^;£^^^). (5.32) 

6. Glassey's argument and blowup in finite time 

Now we are prepared to prove blowup in finite time for negative energy ini- 
tial conditions, by generalizing Glassey's argument familiar from NLS and related 
nonlinear PDE's, to the GP hierarchy. 

The quantity that will be relevant in reproducing Glassey's argument is given 

by 

k 

n(r(t)) := Tr(^xf7('=Ht))- (6.1) 

j=i 

Similarly as in our discussion of the conserved energy, we observe that 

k 

V,{V(t)) = Tr(^x,VHO) 
= fcTr(x?7(i)(i)) 

= kVi{T{t)). (6.2) 

Again, this follows from the fact that 7^'^^ is symmetric in its variables, and from 
the admissibihty of 7^'') {t) for all fc G N, 

Next, we relate d];Vi{t) to the conserved energy per particle. First, let us denote 
by it) the kinetic part of the energy Ei{t) and by E['{t) the potential part of 
the energy Ei{t) i.e. 

E^im) - iTr((-A.)7(i)(i)), 

Enm) = ^TTiBl,_,^^^('+^Ht)). (6.3) 
From (|5.32p . we can relate d1V\(t) to the conserved energy per particle as follows 
dlVxit) = 8 y" duv?^{u;u) + y dX7(X,...,X;X^^^^) 

= \^E^{V{t))^MpEl{V{fy) 

= i6i?i(r(t)) + 4d(p-^)i?r(r(0) 

= 16i?i(r(0)) + 4d/.^:^ /dX7(X,...,X;X,...,X), (6.4) 
P + 2 J , — ■ ■ 
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where we used the fact that Ei{r{t) ) is conserved. 

Now we conclude that for the focusing — —1) GP hierarchy which is either at 
the L-^-critical level (p — Pl'^) or at the L^-supercritical (p > p^a) level, 

d^tViit) < 16^1 (r(0)). (6.5) 

However, the function Vi{t) is nonnegative, so we conclude that if Ei{T{0) ) < 0, 
the solution blows up in finite time. 

To be precise, we infer from (|6.5p that there exists a finite time T* such that 
Vi{t) \ as t / T*. Accordingly, 

1 = Tr(7WW) 

< (Tr(xVi)(t)))V2(Tr( V^t)))^/' 

X 

< C(Tr(a;27(i)(i)))i/2(Tr(-A7(i)(t)))i/2 (6.6) 

where we have first used the Cauchy-Schwarz, and subsequently the Hardy inequal- 
ity. Thus, Tr(-A7(i)(i)) > iVi{t))-^ / oo &s t / T* . 



One can easily verify from (|2.8p that 

k k 

\W'^it)\\nl > ^Tr(~A.^7'''(0) +ETr(-A.,.7('=)(t)) 

= 2fcTr(-A7(i)(i)) y oo (6.7) 

as f y T*. Accordingly, Av^i {T{t)) / cx) as t y T*, which cstabhshes blowup in 
finite time. 
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